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Abstract 

We study QCD-like theories with pseudoreal fermions at finite baryon density. Such 
theories include two-color QCD with quarks in the fundamental representation of the color 
group as well as any-color QCD with quarks in the adjoint color representation. In all such 
theories the lightest baryons are diquarks. At zero chemical potential /i they are, together 
with the pseudoscalar mesons, the Goldstone modes of a spontaneously broken enlarged chiral 
symmetry group. Using symmetry principles, we derive the low-energy effective Lagrangian 
for these particles. We find that a second order phase transition occurs at a value of /i 
equal to half the mass of the Goldstone modes. For values of /i beyond this point the scalar 
diquarks Bose condense and the diquark condensate is nonzero. We calculate the dependence 
of the chiral condensate, the diquark condensate, the baryon charge density, and the masses 
of the diquark and pseudoscalar excitations on /i at finite bare quark mass and scalar diquark 
source. The relevance of our results to lattice QCD calculations and to real three-color QCD 
at finite baryon density is discussed. 

PACS: 11.30.Rd, 12.39.Fe, 12.38.Lg, 71.30.+h 
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QCD; QCD Dirac operator; Lattice QCD; Low-energy effective theory; Chiral Perturbation 
Theory 



1 



1 Introduction 



The study of strong interactions at finite baryon density has a long history. Phenomenological 
knowledge of nuclear forces allows one to obtain a good understanding of equilibrium nuclear 
matter which is relatively dilute. However, these results cannot be simply extended to higher 
densities relevant to neutron stars, supernova explosions and relativistic heavy ion collisions, 
where the microscopic degrees of freedom of QCD, quarks and gluons, become important. 

Understanding QCD at finite density has been a tremendous challenge Unlike finite 

temperature QCD, where significant progress has been achieved using lattice Monte Carlo 
simulations, QCD at finite baryon density is not amenable to such a numerical approach. 
The primary reason is that the determinant of the Euclidean QCD Dirac operator is not real 
at finite baryon chemical potential, /i. 

Recent progress has been achieved analytically by studying QCD at infinite density and 
by studying models of the Nambu-Jona-Lasinio type as well as the instanton liquid model 0. 
It was found that the QCD vacuum at sufficiently high baryon density could become a color 
superconductor 0. In other words, due to the attraction in the color anti-triplet isosinglet 
channel the quark Fermi surface becomes unstable towards the formation of a condensate of 
diquark pairs. 

No first principle lattice calculation methods exist at this moment to study the phe- 
nomenon of color superconductivity. However, the mechanisms which lead to the formation 
of diquark condensates (gluon exchange or instanton induced attraction) can be investigated 
for QCD with two colors. This presents a tremendous advantage. QCD with two colors can 
be studied numerically on the lattice |^ — the determinant of the Dirac operator is real. 
Another class of QCD-like theories with diquark condensation that can be studied on the 
lattice at nonzero chemical potential is QCD with any number of colors and of quarks in 
the adjoint color representation. Clearly, the physics of both types of theories is different 
from that of the three-color QCD, but the differences are easy to understand and classify. 
We hope that understanding the behavior of such theories at finite baryon number density 
will provide us with an additional insight into the phenomenon of diquark condensation in 
three-color QCD. In addition, numerical simulations of both QCD with two colors and QCD 
with adjoint quarks are now being pursued by several groups [^, 

The third class of theories which can be studied in the same way is QCD with the 
phase of the quark determinant quenched. In other words, for each quark such theory 
contains a conjugate quark, with opposite baryon charge. This leads to the appearance of 
colorless diquark states — baryonic pions. The zero-flavor limit Nf ^ of such a theory is 
the quenched approximation of QCD, as has been demonstrated analytically in using a 
random matrix theory at nonzero 

The unifying property of all such theories is the pseudo-reality of the quark representa- 
tions, which manifests itself in the fact that the determinant of the Dirac operator is real. 
Another unifying property (which is ultimately related to this pseudo-reality) is the fact 
that diquarks can make color singlets — these are the baryons of such theories. The Bose- 
Einstein condensation of diquarks, with nonzero baryon charge, may be viewed as baryon 
charge superconductivity (rather than color superconductivity as in three-color QCD). 

In this paper we construct the low-energy chiral Lagrangian describing mesons and 



2 



baryons (diquarks) at finite baryon density for two-color QCD and QCD with adjoint quarks. 
As was pointed out in [§] the /i-dependence of this Lagrangian describing mesons and baryons 
can be fixed by global and certain local fiavor symmetries. Within the domain of validity of 
this Lagrangian we study both its vacuum properties and the mass spectrum of the Gold- 
stone modes. Our goal is twofold: (i) to understand and describe quantitatively the physics 
associated with diquark condensation; (ii) to provide lattice theorists with the qualitative 
and quantitative predictions aiding their data analysis. 

This paper is organized as follows. In section ^ we give an overview of the different 
patterns of chiral symmetry breaking at /i = for three-color QCD, two-color QCD and 
for QCD with adjoint quarks. We relate these patterns to the antiunitary symmetries of 
the Dirac operator, which provide us with a convenient classification in terms of the Dyson 
index j3. Section ^ reviews the global symmetries of the theories, with the emphasis on 
the enlarged SU{2Nf) symmetry. The spontaneous breaking of this symmetry by the chiral 
condensate is the subject of section ^ where we identify the Goldstone excitations. The 
effective Lagrangian at nonzero bare quark mass m is introduced in section ^ In section ^ 
we review the introduction of the chemical potential into the effective Lagrangian by means 
of a local gauge principle p, ^ . The vacuum alignment dictated by the effective Lagrangian 
is analyzed in sections and ^. In section |^ we expand the static part of Lagrangian around 
the minimum and in section |10| we determine the masses of the excitations as functions of fi 
and m from the pole of their respective propagators. In section |ll] we introduce the diquark 
source j and find its effect on the vacuum and the mass spectrum. In section ^ we present 
the dependence of the vacuum condensates and the baryon number density on fi, m and j. 
Finally, in section 13, we rederive the equation of state of the dilute Bose gas of diquarks with 
repulsion and show that it exactly matches the equation of state obtained using our mean 
field analysis of the previous sections. Concluding remarks and discussion are presented in 
section |T3. 



2 Overview and classification 

At zero chemical potential the spontaneous breaking of the chiral symmetry is believed to 
be an essential low-energy property of three color QCD with fundamental fermions. In 
the limit of massless quarks, the QCD Lagrangian is invariant under UiiNf) x UnlNf) 
transformations, but the ground state is not. The analysis of the hadron spectrum and 
numerical simulations on the lattice strongly support this assertion |ll|. The order 
parameter of the spontaneous breaking of chiral symmetry is the chiral quark- ant iquark 
condensate. In the case of two-color QCD with fundamental fermions and in the case of 
any-color QCD with adjoint fermions the symmetry of the Lagrangian is enlarged from 
UilNf) X Uji{Nf) to U{2Nf). In the case of Nc = 2, this symmetry is sometimes referred to 
as the Pauli-Giirsey symmetry [0. Also in these cases there is strong evidence both from 



lattice simulations |T3|, Q and from arguments based on supersymmetry [|T5[ that chiral 
symmetry is broken by a nonzero vacuum expectation value of the chiral condensate. 

The pattern of chiral symmetry breaking is determined by two ideas. The Vafa-Witten 



theorem ||T6| which tells us that vector symmetries cannot be spontaneously broken, and 



the idea of maximum breaking of the axial symmetry [|T^ Il8|. For each of the three classes 
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of theories chiral symmetry is broken spontaneously according to different patterns |18, 17, 
|19|. The axial U(l)^ subgroup of the global symmetry is broken explicitly by the axial 
anomaly in all three cases. The vector-like U{1)b symmetry, corresponding to baryon charge 
conservation, is intact. The remaining symmetries in QCD with three or more colors with 
fundamental fermions are broken according to SUR{Nf) x SUL{Nf) SUv{Nf). For two- 
color QCD with fundamental fermions the symmetry is broken according to SU{2Nf) 
Sp{2Nf), whereas for any-color QCD with adjoint fermions the pattern of symmetry breaking 
is given by SU{2Nf) 0{2Nf). 

We can classify these above three cases by the Dyson index, j3, of the Dirac operator 
with a value of /3 = 2, /3 = 1 and /3 = 4, respectively. The value of (3 is given by the number 
of independent degrees of freedom per matrix element and is determined by the antiunitary 
symmetries of the Dirac operator. It is a concept that originated in Random Matrix Theory 
2l| , ^ , and is important for the Cartan classification of symmetric spaces . 



In the case of two-color QCD the pseudo-real nature of SU (2)coior [|T8|, 0, |T9|, ^ can be 
expressed as the antiunitary symmetry of the Dirac operator T> = juDi, + m, 

[V,T2Cj,K] = or VT2C-f5 = T2C-f5V*, {(3=1) (1) 

where T2 is the color symmetry generator, C is the Dirac charge conjugation matrix and K 
is the complex conjugation operator. Since {T2CKY = 1 it is always possible to find a basis 
in which the Dirac operator becomes real which gives (3=1. 

The symmetry (|I]) persists even at /i 7^ The reality of the Dirac determinant and 



the feasibility of lattice Monte Carlo simulations is the consequence of (|l]). This property 
also allows us to use QCD inequalities at finite fi to show that condensation can only occur 
in the scalar diquark channel [Q. 

In the case of QCD with adjoint quarks the antiunitary symmetry of the Dirac operator 

is 

[V,C-f,K]=0 or VC-f, = C-f,V*. {(3 = 4) (2) 

Since {CK)'^ = —1, it is always possible to find a basis in which the Dirac operator can 
be organized into selfdual (pseudoreal) quaternions. The value of the Dyson index is thus 
/? = 4. 

In both above cases, (3 = 1 and (3 = 4, the antiunitary symmetry leads to enlargement 
of the global symmetry to U{2Nf) (at /i = 0). 

There is no antiunitary symmetry in the case of QCD with three or more colors with 
fundamental quarks. The Dirac operator is a complex matrix, thus (3 = 2. 

In this paper we shall study two classes of theories: (3 = 1 and (3 = 4, aX finite chemical 
potential /i. The main starting point of our analysis is the observation of the fact that in these 
theories the lowest lying baryons belong to the set of Goldstones of the spontaneously broken 
extended flavor symmetry SU{2Nf). As a result the dependence on fi can be described in 



the Chiral Perturbation Theory framework [26, 2/, E 



We shall construct the effective Lagrangian governing the low-momentum modes in the 
theory. These modes are the Goldstone particles of the spontaneously broken global sym- 
metries. The symmetry of the theory is largest at /i = 0, m = 0, and so is the number of 
true Goldstone modes. A nonzero chemical potential /i and/or a bare quark mass m removes 
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part of the symmetry and some of the Goldstone modes acquire masses. Our main goal is 
to find the functional dependence of the masses of such pseudo-Goldstones on fi and m. As 
we proceed we learn about many other properties of the theory: condensates and vacuum 
alignment, phase transitions, etc. 

The two cases: two-color fundamental quarks {(3 = 1) and any-color adjoint quarks 
(/3 = 4) can be analyzed in a similar way. We shall perform this analysis as follows. At each 
step we shall begin with (3=1 case, explaining the concepts and ideas. Then we follow it 
immediately with the same analysis for the (3 = 4 case with emphasis on the comparison 
between the two cases, which will help understand both cases better. We shall use similar 
notations for the objects which are conceptually the same in both cases. In many instances 
we need not rewrite the formulas, only changing the meaning of the notations suffices. All 
the formulas which use any properties specific to either (3 = 1 or (3 = A are explicitly tagged. 
The formulas without such tags are general and apply to both cases. As we shall quickly 
see the two cases naturally complement each other and are, in a certain sense, dual to each 
other. 



3 Global symmetries at m = /i = 
3.1 P=l 

The fermionic part of the QCD Lagrangian with 2 fundamental colors is given by 

We are working in Euclidean space with hermitian 7-matrices and we use spin matrices 
o"^ = {—i,(7k). The quark flavor (as well as color and spin) indices are suppressed and the 
sum over Nf flavors is implied. The symbol D^, denotes color covariant derivative, d^, + iAi,, 
an antihermitian operator, with being a matrix in color algebra, A^, = A'^Ta/2. As 
usual ip = ip^jo = 4'*'^lo and in the Euclidean partition function ip and ip* are independent 
integration variables. 

The fact that the Lagrangian @ has a higher flavor symmetry than the apparent U{Nf) x 
U{Nf) is related to the pseudoreality of the two-color Dirac operator. In particular, the 
conjugate fleld ipn = (J2T24'r transforms similarly (in the same color representation) to ipL- 
Rewriting the Lagrangian (|^) we flnd 

where we have used the well-known properties of the Pauli matrices, — a20"J(T2 = and 
—T2TkT2 = rj, taken into account anticommutativity of Grassman variables, dropped total 
derivatives, and introduced the spinor of dimension 2Nf, 

In this form the U{2Nf) symmetry becomes manifest. Due to the axial anomaly the sym- 
metry in the corresponding quantum theory is only SU(2Nf) (up to discrete symmetries). 
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3.2 /? = 4 

QCD with quarks in the adjoint representation of the color group is described by the La- 
grangian as in but with different notations. The fields ip are now transforming according 
to the adjoint representation of the color group. The color covariant derivative is again given 
by Di, = di, + iAy, but now is given by the antisymmetric matrix {Ai^Y^ = A'^fj^'^, where 
fa^ = /"''^ are the generators of the adjoint representation, i.e., the structure constants. 

The antisymmetric property of the structure constants now replaces the property of the 
fundamental generators: T2TkT2 = — rj = —f^. Using this property we can again 

recast the Lagrangian using spinors of length 2Nf and obtain (^, but with the spinors \E' 
(and ipRj defined by 



:i. - r ■ (-3=4) (6) 





The only difference from (^) is the absence of color T2 matrix in the definition of ipR, and, of 
course, the fact that the spinors ip carry an adjoint, instead of a fundamental, color index. 
Similarly, in terms of the spinors (|^) the SU {2Nf) symmetry of the theory becomes manifest. 

4 Spontaneous symmetry breaking and Goldstones 

4.1 (5=1 

Let us now understand transformation properties of (ipip) 5 the order parameter of the chiral 
symmetry breaking, with respect to the SU {2Nf) symmetry of the theory. We can rewrite 

ipip as 





^^= :n 1 n r =o^'^2r2 , , U/ + h.c.. {(3 = 1) (7) 



The Pauli matrices a2 and T2 ensure antisymmetrization in spin and color indices, to produce 
a spin and color singlet. We see that the condensate is not invariant under all SU{2Nf) 
rotations. The subgroup which leaves (|^) invariant is Sp{2Nf). The Goldstone manifold is 
therefore given by SU{2Nf) / Sp{2Nf) with Nf{2Nf — 1) — 1 independent degrees of freedom. 

Another way of counting the total number of Goldstone modes is starting from the obser- 
vation that the condensate (0) is a product of two fundamental SU {2Nf) flavor representa- 
tions, antisymmetric in flavor indices. Therefore, the condensate belongs to an antisymmetric 
tensor representation the dimension of which is Nf{2Nf — 1). Condensation can occur in 
any of these Nf{2Nf — 1) directions; the fluctuations along the remaining Nf{2Nf — 1) — 1 
directions then become Goldstone modes. 

The effective theory for the Goldstone modes can be written in terms of the fluctuations of 
the orientation of the chiral condensate, E, which, according to the previous paragraph, is an 
antisymmetric unimodular (det S = 1) unitary matrix (exactly Nf{2Nf — 1) — 1 independent 
components). If we denote the equilibrium value of the orientation of the chiral condensate 
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by Sc, the Goldstone manifold given by SU{2Nf) / Sp{2Nf) can be parameterized, according 
to the transformation of S under SU{2Nf), by 



where 

U = exp f ^) and 11 = tt^ ^ — (9) 

The fields vr^ are the Goldstone modes. The implied sum (over a) is over the Nf{2Nf — 1) — 1 
generators of the coset SU {2N f) / Sp{2N f) and in order to simplify the algebra in later 
sections we are using the normalization 

TrX.X, = 2Nf5ab. (10) 

The construction of the Goldstone manifold (^ corresponds to the classification of the 
{2NfY — 1 generators of the SU{2Nf) with respect to a fixed antisymmetric antiunitary 
matrix, in our case given by Sc, into Tk and Xa fl^. The generators leave Ec invariant, 

exp(i0fcTfc)Scexp(z0fcTfc)^ = Sc, i.e. T^T.^ = -T^^T^ . (11) 

By definition, they are the generators of the symplectic group Sp{2Nf). The remaining 
generators, Xa, form the coset SU {2N f) / Sp{2N f) . They obey the relation 

X„Se = SeXj, and thus f/S,f/^ = f/^S,. (12) 

The partition of generators in generators of the coset, X^, and generators of the invariant 
subgroup Sp{2Np), Ti, depends on the matrix T,^- The defining relations ([T^p are left 
unaltered by rotation of Sc according to 

VKV^, (13) 

where V is an SU{2Nf) matrix, and a simultaneous rotation of the generators by 

Xa ^ VXaVl (14) 

This means that the set of broken generators, Xa, changes, if the matrix Sc is changed. 
If we use, as we do below, the following choice for 



-1 

1 



(/3=1), (15) 



the generators Xa can be written in the block representation as 

n=(^'«), (16) 

with Nf X Nf matrices P and Q such that TrP = 0, = P, and = —Q. It is then easy 
to see that the number of independent components in P and Q are, 

Np = Nf-l and Nq = Nf{Nf - 1). (/3 = 1) (17) 
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4.2 /? = 4 



What are the transformation properties of the chiral condensate with respect to the SU{2Nf) 
symmetry in the (3 = 4 case? We can rewrite ipip as 

The Pauh matrix o"2 ensures the antisymmetrization with respect to spin indices, to produce 
a spin singlet. There is no antisymmetrization in color (unlike in the case of fundamental 
colors (3 = 1) — the color singlet is a symmetric product of two adjoint representations. The 
Pauli principle now demands that the SU{2Nf) flavor indices must be symmetric, and they 
are, as is evident in Therefore the condensate belongs to the symmetric (as opposed 

to antisymmetric for (3=1) second rank tensor representation of SU{2Nf), which has 
dimension Nf{2Nf + 1). The condensation can occur in any of the Nf{2Nf + 1) directions. 
The fluctuations in the remaining Nf{2Nf + 1) — 1 directions become Goldstone bosons. 
Alternatively, since the chiral condensate (|TB|) is invariant under S0{2Nf), the Goldstone 
manifold is SU{2Nf)/ S0{2Nf), which gives us the same number of Goldstone modes. 

The Goldstone manifold SU{2Nf) / S0{2Nf) should now be parameterized by symmetric 
unimodular unitary matrices S. The Goldstone fields are introduced in the same way as in 
(1) and (^, i.e., S = UT^cU^, where Sc is now also a symmetric unimodular unitary matrix 
and the implied sum in (y) is over the generators Xa of the coset SU{2Nf)/ S0{2Nf). 

The classification of the generators is also similar to /5 = 1 case. The generators of the 
coset, Xa obey the commutation relation (|1^) with a given symmetric unitary matrix, also 
denoted by Ec, whereas the generators leave Ec invariant as in (0). The Tj are now 
the generators of an S0{2Nf) subgroup (as opposed to Sp{2Nf) in the P = 1 case). The 
remaining generators Xa form a coset SU{2Nf)/ S0{2Nf), which is the Goldstone manifold 
in this case. 

In this case the standard choice for the matrix Ec is 

Sc=(° J). {(3 = 4) (19) 

With this choice, the matrix 11 = rCaXa/ \j2Nf in the coset can be split into 4 blocks of size 
Nf X Nf as in (|TB|). The matrix P is again hermitian and traceless, but the matrix Q is now 
symmetric. Therefore, the counting of the independent degrees of freedom in the (3 = 4 case 
is 

Np = N'j-l and Nq = Nf{Nf + 1). (/3 = 4) (20) 

In both cases, (3=1 and (3 = 4, the kinetic term of the effective Lagrangian describing 
the Goldstone modes should be invariant under the global SU{2Nf) group and under Lorentz 
transformation. The corresponding SU{2Nf) nonlinear sigma-model is given by 

= ^Tr9,Ea,Et, (21) 
where F is the pion decay constant. 
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5 Bare quark mass m 



5.1 [3=1 

If the bare quark mass m (in this paper, the same for all quarks) is not zero, an explicit 
SU{2Nf) breaking term in the effective Lagrangian (|2T1) appears. To determine its form we 
first rewrite the bare mass term using SU{2Nf) notations, 

mijij = ^771^^(72X2 J ) * + ^-c- = ~^''(^2T2M^ + h.c, (/3 = 1) (22) 

where the mass matrix M is given by 

M = mM and M = 1^ J ^ . (/? = 1) (23) 

We see explicitly that the bare mass term is only invariant under an Sp{2Nf) subgroup of 
SU{2Nf). The full SU{2Nf) invariance can be restored if M is also transformed together 
with ^ according to 

^ ^ and M ^ V*MVl (24) 

This extended symmetry must be also manifest in the effective theory, 

S ^ V^V^ and M V*MVl (25) 
The lowest order term induced by the quark mass must therefore have the form 

A.mass = -GReTr(MS) = -777G'ReTr(MS). (26) 

One can view ReTr(ME) as a generalized cosine of the angle between unitary matrices E 
and It is maximal when S is aligned with Therefore the direction of E minimizing 
(p6|) is given by 

Se = M^ (27) 

which leads us to our choice of Eg (pTSl). The mass term comes with a phenomenological 
coefficient, which we denote by G. It is given by the derivative of the vacuum energy with 
respect to m and is, therefore, proportional to the chiral condensate in the chiral limit m 
at /i = (see Section |ir 



G=Mr. (28) 
2Nf ^ ^ 

The resulting Lagrangian with the mass term, 

£^jj = — TrS^Ea^E^ - mG'ReTr(ME), (29) 

is the familiar Chiral Perturbation Theory Lagrangian at lowest order in the momentum 
expansion [E6], |2^. Expanded to second order in the pion fields according to (||), (0), it 
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yields a spectrum with Nf{2Nf — 1) — 1 degenerate (pseudo-)Goldstones with masses given 
by the usual Gell-Mann— Oakes— Renner relation 

= (30) 
We can use this relation to trade G for another parameter, m^, and write 



Tr9^S9^St-2m2ReTr(MS). . (31) 



The symmetry of the theory is reduced from SU{2Nf) to Sp{2Nf) by the mass term. 
Since the chiral condensate does not break any more symmetries in this case we do not have 
any true Goldstones when m ^ 0. 

5.2 /? = 4 

To determine the form of the term in the effective Lagrangian induced by a small bare quark 
mass m we rewrite the quark mass term as 

= lm^^(T2 ~^ ^ ^ + h.c. = -^^^(TsM^ + h.c. , (/? = 4) (32) 

where we have used the spinors of length Nj introduced in section |^. In this case the mass 
matrix is given by 

M = mM with M = ( ° M . (/^ = 4) (33) 



Using the same arguments as in section 5.1 we find that the mass term in the effective 



Lagrangian as dictated by the extended flavor symmetry is given again by (pq), with S now 
being a symmetric matrix and M given by (pBj). Similarly, this term is minimized when 
S = M'^ . With M now taken from (^) we arrive at our choice of Sc (pJ|). 

The symmetry of the theory is reduced from SU{2Nf) to S0{2Nf) by the mass term. 
Since the chiral condensate does not break any more symmetries in this case we do not have 
any truly massless Goldstones when m ^ 0. The masses of all A^/(2A^/+1) (pseudo-) Goldstones 
are equal and are given by the Gell-Mann— Oakes— Renner relation 



6 Chemical potential fi 

In this section we review the introduction of the chemical potential in the effective partition 
function following the approach of This approach relies only on the SU{2Nf) symmetry 
of the theory and the resulting /x-dependent terms are the same for both cases, (3 = 1 and 
(3 = A. 
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6.1 Global symmetries and n 

At nonzero chemical potential the microscopic Lagrangian is given by 

£ = ^jj'yt.D^ijj - ^iijj'yo'ilj + mijjijj. (34) 

As was the case for the mass term, we can also rewrite the baryon charge density in terms 
of the SU {2Nf) spinors. 

(35) 

The physical meaning of +1 and —1 in the baryon charge matrix B is simple: they are the 
baryon charges of the quarks ipL and conjugate quarks ippt- We see that this term is not a 
singlet under SU{2Nf). It transforms in the adjoint representation of this group, in other 
words, the baryon charge is one of the {2NfY — 1 generators of this group. In terms of the 
spinors of length 2Nf, the microscopic Lagrangian is thus given by 



(36) 



where the upper branch corresponds to (3 = 1 and the lower branch to /? = 4. For reasons that 
will become clear in the next subsection, we have introduced the four- vector By = {B, 0). 

As in the case of the quark mass term, the chemical potential term, /i\E'^i?\E' violates 
the SU{2Nf) symmetry. Similarly, we can maintain this symmetry by accompanying the 
rotation of by a corresponding rotation of B, 

^ ^ and 5 ^ VBVK (37) 
Such an extended symmetry must be manifest in the effective Lagrangian, 

S ^ VJ:V^ and B VBVK (38) 
This restricts the lowest order nonderivative term in /i to a linear combination of 

Ij.^Tt{EB^^^B) and ij,^Tt{BB) (39) 

with arbitrary coefficients. Only the first of these terms contains a dependence on the 
Goldstone fields.[] 

At m = the chemical potential breaks the global symmetry of the theory from SU {2Nf) 
down to the usual (as in three-color QCD) SU{Nf)L x SU{Nf)ii x U{1)b- If m 7^ also, 
the symmetry of the theory is only SU{Nf)v x U{1)b- 

^Such type of symmetry breaking terms also occur in the context of the non-hermitian Random Matrix 
Theory @,|29[|o|. 
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6.2 Local symmetry and the coefficient of the fi^ term 



The coefficients of the terms ( pQ]) in the effective Lagrangian can be related to baryon number 
susceptibihty, i.e., the second derivative of the vacuum energy with respect to /i. However, 
unhke G, i.e. the chiral condensate these parameters are not independent. They are 

related to the pion decay constant F by virtue of a local symmetry 

As was observed in p], one can further extend the symmetry (^) of the microscopic 
Lagrangian ( ^61) to include local SU{2Nf) flavor transformations 

^^y^ and ^VB^V^ --Vd^V^. (40) 

l^ 

We can also recover the Lorentz (Euclidean) symmetry by transforming B^ as a four-vector. 

To make such an extended local symmetry (and also Lorentz symmetry) manifest in 
the effective Lagrangian (|3T| ) we must replace the normal derivative by a flavor covariant 
derivative, 

V,S = 9,S-/i(5,S + SfiJ); 

V.St = 5,St + /i(St5, + 5jSt). (41) 

Thus we arrive at the effective Lagrangian of lowest order in the momentum expansion 
F^ 



TrV^SV.St - 2m^ReTr(MS) 
— Tra.S(9.St + 2fiF^TTBJ:^doJ: 



-FVTr (pB^E^B + BB^ - F^m^ReTr (ifs) . (42) 

It is important to note that the dependence on /i comes with no additional parameters. It 
is completely fixed, by the local symmetry, in terms of an already existing parameter F. 



7 Vacuum alignment 



The static part of the Lagrangian (^2]) determines the vacuum alignment of the field E as 
well as the masses of the excitations. This part of the Lagrangian has the form 



Cst{^) = -FVTr (JIB^E^B + BB) - F2m^ReTr(MS) 
F'^ml 



2 

— Tr (t.B'^T^B + BB) - 2ReTr(MS) 



(43) 



where we introduced x = 2/i/m^. The and terms in (|43|) compete for the direction 
of the condensation which we denote by S. For a; = the orientation of S is determined by 
the mass matrix S = . This value we denoted by S^. This is the orientation of the usual 
chiral condensate which carries no baryon charge. When a; = oo p[ the static Lagrangian is 
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minimized on a manifold in the space of S, from which we choose the following value and 
denote it by S^, 

(/? = 1) (44) 

where / is an antisymmetric Nf x Nf matrix, which, written in {Nf/2) x {Nf/2) blocks, has 
the form 

(45) 

This minimum (even for m ^ 0, unlike Sc) has a degeneracy, which for Nf = 2, is simply a 
rotation by the generator B. The condensate breaks spontaneously the baryon number 
symmetry. The degeneracy leads to the appearance of a massless Goldstone boson. For 
Nf > 2 this condensate breaks more than just the U{1)b symmetry. For m = it also 
breaks SU^Nf)^ x SU{Nf)ji down to Sp{Nf)i x Sp{Nf)ji given by the rotations leaving the 
block matrix / in ( ^if ) invariant p|. When m 7^ it breaks the SU{Nf)v symmetry down 
to Sp{Nf)v.f\ 

The situation is similar in the (5 = 4 case. The only difference is that now the matrix S 
is a symmetric SU {2Nf) matrix. It rotates between the value of Sc given by (0) and S^, 
which we chose as 

{(5 = 4) (46) 

The discussion of the preceding paragraph carries over to the (3 = 4 case with the already 
familiar substitution of Sp{Nf) by SO{Nf). The condensate breaks SU{Nf)L^ji flavor 
symmetries down to SO{Nf)L ji, in addition to breaking U{1)b- 

At intermediate values of x the orientation of the condensate S rotates, as a function of 
X, from Sc to T^^i. We shall prove that it can always be written as the linear combination 

S = = Sc COS a + Srf sin a, (47) 

where a = for x = and 7r/2 for x = 00. 

The angle a is a function of x. It is determined by substituting the value of S given by 




(^) into (^31) . This results in the static Lagrangian 

x'^ 

— (cos 2a — 1) — 2 cos a 



At(S,) = F^mlNf 



Minimizing with respect to a we find 
a = 0, when x < 1; 



(4J 



cos a = — , when x > 1. (49) 



X^ 



^Similarly to the alignment of Ec to the bare quark mass matrix . the direction of is determined 
by an external diquark source, J, which, as we shall see in section |ll|, breaks the degeneracy. When J is 
zero, we can choose any orientation of E^ within the manifold of minima, and the results will not change 
due to the symmetry relating all such minima. 
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Note that a = is always an extremum of Pq) but it becomes a maximum for x > 1. 

We find that the condensate is a non-analytic function of x = 2/i/m^. There is a second 
order phase transition as a function of /i at x = 1. For small /z < the vacuum does not 

change. At n = a transition occurs and the direction of the condensate starts rotating 

from that of Sc to that of S^. A nonzero value of the projection on E^, proportional to sin a, 
means that the vacuum breaks the baryon number symmetry spontaneously. This happens 
at the value of equal to 1/2 of the mass of the lightest baryon in the system - the di quark, 
or, the baryonic pion. 

8 Global minimum 

In this section we show that the minimum of the static potential is given by defined 
in (0). The argument in this section shows that it is an absolute, or global, minimum 
of £st(S). In the next section we expand the Goldstone fields to second order about this 
minimum. 

8.1 P=l 

We decompose the antisymmetric unitary matrix S into 4 blocks of size Nf x Nf, 

{(3 = 1) (50) 

The antisymmetric matrices A and B satisfy the following unitarity constraints: 

AA^ + CC^ = 1, BB^ + {C^Cf = 1, AC* = CBK {p = 1) (51) 

Taking this into account we can express Cst entirely in terms of the matrix C, 




(52) 



Ignore, for the moment, the constraints on the matrix elements of C. The trace in 
can be viewed as the distance in the 2Nj dimensional space of real and imaginary parts of the 
matrix elements of C from the point given by diagonal matrix 1/x^. For x > 1 the absolute 
minimum is achieved when C = 1/x^. The matrices A and B can be chosen, for example, 
as ^ = -B = ilJl — 1/a;^, to satisfy all the constraints (pTI). If we define cos a = 1/x^, we 



observe that the resulting matrix in (|3) is given by Sq (|47|) . 

When X < 1, we notice that unitarity constraints (|51|) demand that TrCC^ < 1. This 
means that we have to consider the points in the space of C only within the unit hypersphere 
around C = 0. It is easy to see that the minimum distance within this sphere is at the surface 
point closest to 1/x^, i.e. at C = 1. The constraints (0) can be satisfied only hj A = B = 



and the resulting matrix in (|50|) is Sc, the minimum of the static potential at fi = 0. 
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8.2 /? = 4 



The same analysis apphes in the (3 = A case. Since S is now a symmetric matrix its block 
decomposition is given by 



A +C 
B 



{(3 = 4) (53) 



with symmetric matrices A and B and with unitarity constraints 

AX^ + CC^ = 1, BB^ + {C^Cf = 1, AC* = -CBl {(3 = 4) (54) 

The static part of the Lagrangian is the same as in (0). The minimum is given by C = 1 
for X < 1 and C = for x > 1. The matrices A and B are not determined uniquely by 



the unitarity constraints. One can check that A = B = i^l — (as in (|47|) ) satisfy these 
constraints. 



9 Curvatures at the minimum 



We have established that the global minimum of the effective Lagrangian is achieved at the 
value of the field S given by (j^7\} and (^). The orientation of the condensate rotates as a 
function of x = 2fi/m.,^ (nonanalytic at x = 1). In this section we shall expand the effective 
Lagrangian up to the second order in the fluctuations of S which will help us to determine 
the (pseudo-)Goldstone masses and their dependence on x. 



9.1 Normal phase: fi < 

When X < 1 the vacuum orientation of S does not depend on x and is given by S = Sc. 
Expanding S around Sc using the Goldstone fields defined in (P) according to 

j: = UKU^ = u'K=(i + f-^, + ..)K, (55) 



As in the case of /3 = 1, the case of odd number of flavors has to be investigated separately. Let us 
consider Nf — 1. Then a symmetric unitary matrix has two degrees of freedom and can be parameterized as 



ismUe^ cost 
cos 6* i sin 6*6 



Notice that the determinant of E has to be equal to the determinant of the mass matrix M (|3^). We thus 
have that the matrix C in ( ^3|) is just a number parameterized according to C = cos 9, leading to the effective 
potential (up to constants) 

x^Tt (c^^X) (c-X) - x^ - \ ^ x^ cos"^ 9 - 2cos9 - x"^. 
\ x-' J \ x^ ) x^ 

The minimum is at = for x < \ and at cos0 ~ l/x"^ for x > 1. In the diquark condensation phase 
{x > 1), we find one massless excitation (the 0-mode), a nonzero baryon density and nonzero chiral and 
diquark condensates. 
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we find 



X 



^Tr[5,n]2 + Trn2 



+ . . . 



(56) 



where ellipsis denotes terms of higher powers of 11. The commutator in this Lagrangian can 
be written as 



[B,U] = bU, 



(57) 



with b the baryon charge of the pseudo-Goldstone modes. Since all our pseudo-Goldstone 
modes are quark-antiquark or (anti-)diquark states the values of b are 6 = 0, ±2. The 
curvature of the P, Q and ([T6| ) modes, thus depends on through the baryon charge and 
is given by — (fo/x)^. For example, in the /3 = 1 case of Nj = 2, there are 3 pseudoscalar 
mesons, 1 diquark and 1 antidiquark. Using the block decomposition of the generators (|16D 
we find 



/:st(S) = /:st(Se) + m^Tr P2 + (1 _ ^2^QQ 



it 



+ . . . . 



(58) 



Two comments are in order here. First, note that, the curvature (and the mass) of the 
meson modes P does not depend on fi in the normal phase. Second, at x = 1 the curvature 
of the diquark modes Q vanish, which signals a phase transition and the onset of diquark 
condensation. 



9.2 Diquark condensation phase: fi > 

In this phase the condensate begins to rotate according to (^ 
also written as 



T 



V^J:,, where 



(^9]). This rotation can be 

(59) 



and X2 is the generator that rotates Sc into S^. Comparing 
Erf = iX2^c- 



we find 



(60) 



This generator belongs to the set of broken generators with respect to (as well as with 
respect to E^) since it satisfies (|T2|). 

We could parameterize fluctuations of S around the vacuum value given by as 



T 
a 1 



(61) 



where Ua are unitary matrices generated by rotated (in the sense of (p^)) generators, V^X^V^, 
instead of Xa (|^), (§). However, as we shall find, the meson mass matrix is diagonal in the 
basis given by the parametrization 



(62) 



where we have rotated the coset generators Xa back to their a = values using (|T^, ( P^^D 
and (|5|). 
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Before expanding the static Lagrangian, we substitute (^) into (^3[) to obtain 



yTr f f/^S, (vJBVo)" j:1 (u^Y (vJBVo) + BB 



T 



-2ReTr (iv^ MV^) f/^S 



(63) 



The rotated values of B and Sc (by angle —a in the sense of ( ^Of ) ) can be expressed as 



VlBV^ = V.aBVl^ = Bcosa + Xi si 



sm a 



V^MVa = (VI^MVT^)^ = SL„ = Sjcosa - S^sina, 
where we have introduced Xi by 
Xi = iBX2. 



(64) 



(65) 



This generator, similarly to X2, belongs to the coset of broken generators with respect to 
(0). It is the generator into which B is rotated while is rotated into E^^. 
Substituting (Q) into (^) and expanding to second order in 11 (|55| ) we find 



At, (SI 



X 



sin 2a + sin a 

2 



Tr {X2U) 



+ - 



mi 



^ (Tt[B, cos^ a - Tr[Xi, n]^ sin^ a) + TrH^ cos a 



+ 



(66) 



As should be expected, the linear term vanishes due to ( ^9] ) and we shall concentrate now 
on the quadratic term. 



9.2.1 p = l 

At this point we need the explicit form of Xi for (3 = 1. According to (|60|) we have 

X2 = -t^d^l = 
For Xi we therefore find 
Xi = iBX2 = 



/ 
-/ 



il 
il 



(/5=1) 



(67) 



(68) 



We then use the block decomposition of 11 ([iq) and find 

£st(S) = £st(S„) + m^Tr [-x^ [QQ^ cos^ « - (P| + gflQj,) sin^ a) 



+ {P^ + QQ^) cos a 
We have introduced the following projections of P and Q: 



+ ... 



iP = 1) (69) 



Ps = ^{P + IP'' I) and = 1 (P - JP^J) ; 

Qfl = ^(Q + /Q^/) and g, = i(g-/gt/). 



(/5 = 1) 



(70) 
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These projections are orthogonal, TtPsPa = R-eTrQ/jQ| = 0. Using this fact and the 
relation cos a = we obtain 

£st(S) =£st(Sa)+m^x2Tr [Q^Qj^sin^a + Plcos^a + Pj] +.... {P = 1) (71) 

From (|7TD we can now read off the curvatures for the different multiplets of the (pseudo- 
)Goldstone modes. We see that there is a true fiat direction, Qj, which describes true 
massless Goldstones of the diquark condensation phase. These fields are the phases (a single 
f/(l)ij phase for Nf = 2 and a set of f/(l)B x SU{Nf)v/ Sp{Nf)v phases for other Nf) of the 
diquark condensate. As we shall see in the next section, the linear derivative terms in the 
effective Lagrangian mix Qj and Qn, thus the actual true Goldstone excitations are certain 
linear combinations of Qj and Qr. 

The degeneracies of the multiplets follow from the definitions ( [70D and are given by 

No, - N^.^^^I^. (P^l) (72) 

They correspond to representations of the group Sp{Nf), which is the residual symmetry 
remaining intact after spontaneous breaking. The number of the fiat directions Qi exactly 
matches the number of the broken generators in SU{Nf)v x U{1)b Sp{Nf)v- 

9.2.2 /? = 4 

The derivation of the curvatures of the effective potential Cgt in the /? = 4 case follows the 
same lines as in the = 1 case. All the differences stem from the fact that S is a symmetric 
unitary matrix in this case. Consequently, Sc and are given by (|T9|) and (^) instead of 



(p!5|) and (^if ). In particular, the generator which rotates Ec into is now given by 

X2 = -^S,St = f ° i V (/3 = 4) (73) 



and the generator Xi, into which B rotates and which is responsible for the splitting of both 
P and Q branches (in the block decomposition (0)), is given by 



X, = tBX,= ^_^ oj- (/3 = 4) (74) 

It is now straightforward to rewrite the quadratic part of the Lagrangian (pB|) in terms of 



the PQ-block decomposition. The mass matrix becomes diagonal in terms of the following 
projections of P and Q: 

Ps = ^{P + P^) and Pa = 1{P-P1; 

Qr=\{Q + Q^) and Q, = i(g-gt). (^ = 4) (75) 
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The static part of the effective Lagrangian in terms of these fields is given by 

At(S) = £st(S„) + mlx^Ti [QnQi sin" a + PI cos^ a + Pj] + . . . , (/3 = 4) (76) 

which is similar to (0), but S and A labels exchange places. 

The main difference is in the dimension of the degenerate multiplets which matches the 
representations of the corresponding residual symmetry group. By inspection one finds 

No, - N„^Mk±ll, (^.4) (77) 

10 Spectrum 

In order to complete our program and determine the spectrum of low-lying excitations we 



must take into account the derivative terms in the Lagrangian (^). They contribute in 
a non-trivial way to the quadratic form in the Goldstone fields. Due to the non-Lorentz 
invariant nature of the system we study at finite fi, the dispersion laws0 do not have a simple 
form, E"^ = + m?. The mass, as measured on the lattice and given by the exponential 
fall-off of the propagator at large Euclidean time, is the value of E, i.e. ipo? at the pole of 
the propagator at p= 0. These pole masses, or rest energies, we shall now evaluate. 



10.1 Normal phase: /i < 

Expanding the derivative terms in the Lagrangian (^2|) in the same way as we expanded the 
static part in the previous section we obtain 

C{T.) = + ^Tr {(9,n - n])' + mlW] + . . . , (78) 



where By = BS^q, as defined earlier in Section |6.2| . Using the Fourier decomposition of 11 



n(a;) = Ilpexp{—ipx) = lip exp{—Exo — ipx), (79) 
p p 

we find that the dispersion law has the generic form 

{E + bfif = p^ + ml, or E= -bfi + ^p^ + ml, (80) 

where h is the baryon charge of the given excitation. This is in agreement with the fact that 
the effect of /i on each state is simply an energy shift by — In particular, the dispersion 
law of the P-type Goldstone modes, which carry no baryon charge, are not affected by yU at 

"^In Euclidean field theory the dispersion relation is given by the poles of the propagator in the E = ipo- 
plane. 
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all. The rest energies of the diquarks Q and antidiquarks are shifted according to their 
charges. To summarize, the dispersion laws in the normal phase are given by 

P : E=(ml + p')"^ ■ 

Q : E = {ml + p^y^^ - 2fi. (81) 

The masses of these excitations as well as their degeneracies are given in Tables |I] and ^ 
A schematic picture of the mass dependence is shown in Figures |^ and ||. 



10.2 Diquark condensation phase: fi > 

In this phase the ground state changes and we need to expand around the rotated value of 
the condensate E = Sq. We find^ 

£(S) = /:(S„) + ^TT{[d,U-fi[V^B,VlU])[d,U-fi[V^B,V^,U]) 

+ mlU^ cos a} + .... (82) 
Expanding the product we find 

C{E) = ^Tr {(9,n)2 - 2ficosa[B,U]doU} + C^U^) + ... (83) 

with the static part given by (|66[). We observe, first of all, that the linear derivative term 
contains only the charged fields, i.e. Q and Q^. The dispersion laws of the P fields remain 
unaffected by the linear term. It remains Lorentz invariant, with mass given by the curvature 
of the static part of the Lagrangian. In order to determine the dispersion laws for the Q and 
fields we need to solve a secular equation obtained by substituting Fourier decomposition 
of Q's into 

/:(E) = C{E^) + TT{[d,Q^jid,QR + d,Q\d,Qi) 

-Afi cos a [Q\doQR + QidoQi) + VqJ^Q^ sin^ «} + ..., (84) 

where we only wrote the Q-dependent terms. This expression is the same for either (3=1 
or /3 = 4; only the definition of Q's is different: eqs. (|70D or eqs. ([75|). The secular equation 
has the form 

^ / E"^ — —AfiE cos a 

I —4:fiE cos a E"^ — p^ — 4/i^ sin^ a 

^Note that though the mass of the diquarks Q vanishes at the transition point, /i = 171^^/2, the dispersion 
law is not linear, but quadratic, E ^ p^. This is related to the well-known critical slowing down at a second 
order phase transition. 

useful observation which helps to write this in such a compact form is that Tr[Xi,n]9on = for any 
n. This can be checked explicitly using the block decomposition of 11, but it is easier to see that once you 
realize that the commutator of two X-like generators is a T-like generator, and its trace with another X-like 
generator, such as is zero. 
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We see that Qr and Qi are mixed by the hnear derivative term. The mixing is maximal 
when a = 0, and vanishes when a = tt/2, i.e. when fi ^ 171^^. For any a there is a solution 
for which the rest energy -E'(O) = — the Goldstone boson branch, which we denote Q. For 
a = it is given hj Q = Q = Qr + Qi (at rest, p = 0). For a = tt/2 it is entirely Qi. The 
other solution of the secular equation, the linear combination of Qr and Qi orthogonal to 
Q, is massive. We denote it by 

The dispersion laws for the meson multiplets Ps, Pa and the mixed diquark-antidiquark 
multiplets Q and are given by (in the (3 = 1 case) the following relations: 

Ps : = p2 + mlx^; 

D 7712 „2 I 2 2 2 

Pa ■ E = p + m^x cos a; 

: E'^ = p^ + 2i?{l + 3 cos^ a) + 2/iy'/i2(l + 3 cos^ af + 4p2 cos^ a; 
Q : ^2 = p2 + 2/i2(l + 3 cos^ a) - 2^1^ + 3 cos2 a)2 + 4p2 cos2 a. (/5 = 1) (86) 

The only difference in the dispersion laws for /? = 4 is the interchange of Ps and Pa-, and 
also that the degeneracy of the multiplets as given by (|72D and is different. Note that 
the dispersion relation for the lowest lying mode, i.e. Q, is linear: E p. The slope is a 
function of /i: it vanishes at the transition point, x = 1, and approaches 1 for large /i. The 
linear slope of low energy excitations is characteristic of superfluidity. 

The pole masses of the excitations, i.e. the position of the pole at zero momentum,£'(0), 
are given in Table |I| for (3 = 1 and Table ^ for (3 = 4. The representations of the residual 
symmetry groups are denoted by their Young diagrams. Because mesons are quark-quark 
or quark-antiquark pairs they transform as rank two tensors (unless they are singlets). The 
representations can be uniquely identified by their dimensions found in (|1^, ( [72D and (pO|), 
(frzl). An explicit form of the representations is given in the Appendix. Figures and ^ show 
schematic pictures these results, together with the residual symmetry groups in the different 
phases. In particular, one observes that the spectrum is continuous at the transition point 
X = 1. 

11 Diquark source 

In three-color QCD the diquark condensate is not a color- singlet, therefore, one cannot study 
the phenomenon of diquark condensation by applying an external diquark source: such a 
source term will not be gauge invariant. In the theories which we study in this paper, e.g., in 
two-color QCD, the diquark condensate is colorless. Therefore one can add a gauge-invariant 
source j to the theory. This source plays a role similar to the role the quark mass m plays 
with respect to the chiral condensate. Such a non-zero source term j is in fact used in lattice 
simulations in order to measure the value of the diquark condensate in the limit j — >■ 0. 
In this section we show how the results of the previous sections are modified in the presence 
of a non-zero j. 
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Sp(2Nf) 



SU(Nf) X U(l) 



Sp(Nf) 



Sp(Nf) X Sp(Nf) 

(□,□) 



m 



P: \*i 




Q: 1 + 



(1,1) + (1,H) + (H,1) 



x=- 



Figure 1: Spectrum of two-color QCD (/5 = 1) at finite and m (schematic). The branches 
are labeled according to (|8T|), ( |70| ) and ^) and by Young's diagrams showing the multiplet 
structure under the corresponding residual symmetry. The star denotes the conjugate fun- 
damental representation. The residual symmetry groups are marked above the plot. The 
degeneracies are given by (|17D and ([72|) . For example, for the case of Nf = 2 we find Np = 3, 
Nq = Ngt = 1, and Np^ = 3. The branch P4 does not exist for Nf = 2. On the left /right of 
the plot the residual symmetries and multiplet structure of the corresponding limiting cases 
fi = and m = are shown. 
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Table 1: Spectrum (masses and degeneracies) of two-color QCD with fundamental quarks 
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Table 2: Spectrum (masses and degeneracies) of any-color QCD with adjoint quarks (/3 = 4). 
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SO(2Nf) 



SO(Nf) 




SO(Nf) X SO(Nf) 

(□:□) 



(i,i) + (i,m) + ([ 



1,1) 



Figure 2: Spectrum of any-co\oT QCD with adjoint quarks {(3 = 4) at finite fi and m 
(schematic). The branches are labeled according to (|8TD , ( [70l) and (|75D and by Young's 
diagrams showing the multiplet structure under the corresponding residual symmetry. The 
residual symmetry groups are marked above the plot. The degeneracies are given in Table 
|. For example, for the case of Nf = 2: Np = 3, Nq = 6, Nq = Nq^ = 3, Np^ = 2 and 
Np^ = 1. On the left /right of the plot the residual symmetries and multiplet structure of 
the corresponding limiting cases n = and m = are shown. 
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11.1 [3 = 1 



Let us first rewrite the scalar diquark source term in the macroscopic theory using our 
SU{2Nf) spinor notations. In the notations of section ^ we find 

- i|v^C75r2/V' + h.c. = -|^^fT2r2 |^ ,7 ) * + h-^' = -^^^^^2^2^^ + h.c, 

J = jJ and i=('d °/)' = ^^^^ 

the antisymmetric matrix / is defined in (^51) , and the summation over the Nf flavor indices 
has been suppressed. Comparing with the quark mass term (p^ we find that the two belong 
to the same multiplet (i.e. they transform into one another) under the SU{2Nf) rotations.]] 
We can write the sum of the diquark source and the bare mass term as 

m^^ - i^-{^^Cj,T2lij + h.c.) = -^^^a2r2M^^, (/3 = 1) (88) 

where 



= mM + jj = + j2(Mcos0 + Jsin0) = ^Jm'^+fM^. (89) 
We have used a notation similar to { ^7\) with 



tan0=— . (90) 
m 

It is easy to see how the diquark source term modifies the effective Lagrangian (^2]): we 
need to replace the mass matrix rriGM in ( pPj ) by ^/m'^ + j'^GM^. The Gell-Mann— Oakes— Renner 
relation at /i = becomes 

< = (91) 

Using this fact we can write the effective Lagrangian for the theory with both mass term 
and diquark source (compare to (0)) as 

r - ^ 

CesiX) = — TrV^SV^S^" - 2m^ReTr(M<^S) . (92) 



Now we can repeat the steps we performed in section H, but for the Lagrangian (|9^) . 

Since the source term we introduced favors the direction of S^, we expect that the 
minimum of the Lagrangian (|9^) , S, is again given by a linear combination of and as 
ini3), 

S = E„ = cos a + Ed sin a = VJlcV^, where K = e*"^^ , (93) 



^There is a freedom in the choice of the orientation of the diquark source J. It is precisely the same 
freedom as the one corresponding to the SU{Nf) x U{1)/ Sp{Nf) degeneracy of the diquark condensation 
vacuum S^. An infinitesimal diquark source j determines the vacuum orientation of E^. Our choices reflect 
this fact. Indeed = J^. 
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and X2 is the generator that rotates into as before. The value of a, as determined 
by the saddle-point equations, now depends on the value of the diquark source. The proof 
of Section |^ is not easily extendable to this case. We shall continue on the assumption that 
the minimum is global. We shall prove, however, that it is a minimum when we expand to 
second order in fluctuations of E and find no linear terms and a positive quadratic form. If 
we substitute this ansatz into the effective Lagrangian (|92D we find 



F^mlN 



X 



— (cos 2a — 1) — 2 cos(a 



(94) 



The dependence of the angle a on x follows from minimizing this Lagrangian and is given 
by 



X cos a sm a 



sm a 



(95) 



which is different from (|^), but coincides with it in the limit — 0, of course. For finite 
the angle a is already non-zero (a = 0) at x = as it should be, since a diquark source drives 
a non-zero diquark condensate. At x = 00, the value of a is 7r/2 independently of 0. The 
major difference from (^) is that the dependence of a on x is analytic. There is no phase 
transition as a function of /i when the diquark source j is nonzero. This is to be expected: 
the diquark source plays the role an external magnetic field plays in a ferromagnet. The 
diquark source explicitly breaks SU{Nf)v x U{1)b symmetry down to Sp{Nf)Y and there is 
only one phase (the one with residual Sp{Nf)v symmetry) for all values of x. 

Our next step is to expand in powers of fluctuations of S around S. This can be done 
following the algebra in section ^ with minimal modifications. The Goldstone manifold is 
again parameterized by 



but the second identity in 



is now given by 



V;^M^Va = (V-aMlV^^y = = Nicosia - 0) - S^sinfa - 



(96) 



(97) 



whereas the first identity is the same as before. Using this algebra, the effective Lagrangian 
to second order in the pion fields can be expressed as 



X cosasma — sm a 



/:(s) = c,,ij:^) + Fml 

+ -Ti[{d^Uf -2ficosa[B,U]doU} 



Tr(X2n) 



+ 



mz 



X 

T 



(Tt[B, cos^ a - Tr[Xi, H]^ sin^ a) + TrH^ cos(a 



+ . 



(98) 



Inspecting the linear terms we find that they vanish if a is chosen according to (p5|). This 
means that our ansatz (p3|) is indeed an extremum. 

In order to find the mass spectrum we have to study the term of second order in the 
H-fields in the Lagrangian (^), and determine the dispersion laws of the different (pseudo- 
) Goldstone fields. Using block decomposition of the generators 11 ([T6|), projections ( ffOD of 
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P and Q and the relationship (|95|) between x and a we find 



£st(Sa) + Tr {d,Q^j,d,QR + 9,g}9,g/ - 4/icosa g}9ogi? + g}j9og/ 



g/gl^-^ + g^gi? sm^ « + - 

sm a ' 



sm I 



sm a 



+Tr 



O^PaO^Pa + Pim^ x' cos' a + - 



sm - 



sm (p 
sin a 



sm a 



+ 



(/5 = 1) 



(99) 



As in the case of a vanishing diquark source studied in the previous Section, the hnear 
derivative term contains only the charged fields Q and g^ Therefore the dispersion relations 
for the neutral fields P are not affected by this term: the masses of the P fields are given 
by the curvature of the static part of the Lagrangian at the minimum. However the linear 
derivative term mixes the charged Goldstones. In order to obtain the dispersion laws for the 
g fields, as in Section 10, we have to solve a secular equation obtained by substituting the 
Fourier decomposition of these fields into the effective Lagrangian (|99|) . The secular equation 
reads 



det 



m 



2 sin 4> 



P """"sin a 

—4:fiE COS a 



—AjjiE cos a 



_^2sin^ 



0. 



(100) 



T sin a 



The dispersion laws for the different Goldstone modes are therefore found to be given by 

^ sin d) \ 



Ps 
Pa 



E^ 
E' 
E^ 



p [x + 



sm a 



p + X cos a + 



sm I 



sm a 



p 



sm( 
sin a 



+ 2/2^(1 + 3 cos^ a) 



g 



E' 



+2/iy /i2(l + 3 cos^ a)^ + 4 cos^ a {p^ 
sin(' 



sm I 



mt 



sm a 



9 2 

p" + m. 



sm a 



+ 2/2^(1 + 3 cos^ a) 



sm I 



-2/iW/i2(l + 3cos2a)2 + 4cos2a (p2 + ^2__i). (/5 = i), (lOl) 



sm a 



where a is an implicit function of x and (j) given by (p5|). Note that the Q are no longer 
true Goldstone modes, since the symmetry, which is broken spontaneously in the diquark 
condensation phase x > 1 at j = 0, is now broken explicitly by the diquark source term. The 
masses of the different multiplets are given by the value of at p = 0. The positivity of all 
masses shows that the minimum given by the saddle-point equation ( PBj ) is at least a local 
minimum. This mass spectrum at a small non-zero diquark source is depicted in Fig. R. 
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11.2 (3 = 4 

In the (3 = 4 case the diquark source term has the form^ 

j.^T / i \ - . 1 



i-i(^ C'y^ip + h.c. = -^^(T2 



J = jJ 



and 



J 



. ^ 

i 
i 



\I/ + h.c. 



^^^(TaJ^ + h.c., 



(/5 = 4) 



(102) 



The analysis follows the same lines as in the (3 = 1 case. We again introduce the combined 
mass matrix 



= mM + jJ= Jrn^ + pM, 



(103) 



with M and J defined for /3 = 4 in (|33D and ( p,02|) . In the expansion of the Lagrangian to 
second order in the Il-fields for /3 = 1 we have only used the general commutation properties 
of M(j, which, with a proper redefinition of the generators Xi and X2, are the same in the case 
(3 = 4. We thus obtain the same second order Lagrangian (|98D as for (3 = 1 resulting in the 
same saddle point equation ( P5| ) and the same mean field value of the free energy. However, 
the mass spectrum is slightly different. We express the quadratic part of the Lagrangian in 
terms of the block-decomposition (|^) and the projections (|75|). The final result is 



CJ^o,) + Tr 



(d^Q^^d^QR + d^Q\d^Qi) - 4/icosa {Q\doQR + QrOoQi 



-l-m^Tr 



QiQi—^ + QrQr sm^ a + - 
sm a 



sm I 



sm a 



+Tr 

+Tr 



d^Psd^Ps + Plml cos^ a + 



sm I 



sm a 



d.PAduPA + Plml \x' + - 



sm I 



sm a 



+ .... 



(/? = 4) 



(104) 



As in the (3=1 case, the linear derivative term mixes the Q and fields. In order 
to get the dispersion relations of these fields, the same secular equation as before ( |100|) has 
to be solved. The dispersion laws are the same as in ( |101[ ) with the familiar exchange of 
labels S ^ A, and the degeneracy of the multiplets given by (^) instead of (|72D . The 
corresponding mass spectrum at a small non-zero j is shown in Fig. 13. 



12 Condensate and baryon charge densities 

In this section we derive the dependence of the chiral condensate, the diquark condensate 
and the baryon charge densities. These densities are obtained by differentiating the vacuum 
energy density with respect to m, j and /i, 

m = -^-^; m = -^-^; ns = -^. (105) 



*As before, our previous choice of corresponds to our choice of J, i.e., = J^. 
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Figure 3: Mass spectrum for a small non- vanishing source j = 0.1m. 
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1 case and of 



In the effective theory the vacuum energy density is given by the value of the effective 
Lagrangian (^21) at the minimum, 



The shorthand symbol ijjip denotes the l.h.s. of the equation (87) in the /? 
( |102|) in the /? = 4 case without factor j. 



-2fi^TT 5 + BB) - 2m^ReTr (e^E„) 

= —ANfF'^ii^ sin^ a — 2NfG{m cos a + j sin a) , 
where we used (|90|) and (^T]). Differentiating]^, we find 

{ipip) = 2NfGcosa; (tpij) = 2NfGsma] ub = SNfF'^fisin'^ a, 



(106) 



(107) 



where the angle of the rotation of the condensate, a, depends implicitly on /i, m and j 
through the solution of (pS]). Note, in particular, that the sum {ipip)^ + {ip'ip)^ does not 
depend on fi, m or j (to the order in Chiral Perturbation Theory we are working), which 
reflects the fact that the condensate rotates. 

For j = the dependence of a on /x is simple: a = 0, for yU < m^/2, and a = 
arccos[m^/(4/i^)] otherwise (^Sp. That means the densities ( |107| ) are constant for /i < m^/2, 
as they should since the vacuum state does not change until /i reaches mT,/2. The results 
for the condensate and baryon charge densities are summarized in Table 3. 



Table 3: The values of the chiral condensate, {ipip), the diquark condensate, {ipip), and the 
baryon density in the two phases of the theory. 



phase 








/i < m^/2 
/i > m^/2 


(#)o 







S^,NfF\l-{^)') 



For j 7^ the angle a varies between the values and 7r/2 as a function of ^ and m, 
determined by the saddle point equation (p5|) a;^ sin a cos a = sin(a — (p). The saddle point 
is a quartic equation which, in principle, can be solved analytically for arbitrary j. The 
dependence of the densities ( |107|) on ^ is shown in Fig. ^ for j = and in Fig. ^ for j ^ 0. 



13 Equation of state of the non-ideal diquark gas 

We now take a different look at the system we study (quarks with two colors or any-color 
adjoint quarks) near the threshold, /i ~ m^/2. First of all, we observe, that for ^ very close 
to m7r/2 the density of baryons is very (arbitrarily) small. This means that we can describe 
the system by an almost ideal gas. It should be a Bose gas, since the baryons, i.e. diquarks, 
are bosons. The equation of state of a dilute non-ideal Bose gas is a textbook problem 

^The implicit dependence of a on m, j and is not contributing because dE^a-c/da = 0. 
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Figure 4: The magnitudes of the chiral (V'V') and the diquark {tp'ip) condensates in units of 
{'ipip)o — 2NfG as a function of ii/uIt^ for zero diquark source. Also the density of the baryon 
charge in units of S2NfF'^m^ is shown. 
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Figure 5: The magnitudes of the chiral (V"^) and the diquark (^^) condensates in units 
of 2Nf{'^'ip)o as a function of ///m^ at small non-zero diquark source j — 0.1m. Also the 
density of the baryon charge in units of 32A^jF^m^ is shown. 
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5T|]. The interaction between the bosons is crucial here (the dependence on its strength is 
non-analytic, e.g. the attractive Bose gas is unstable). The baryons in our system repel 
each other, as we shall see. Moreover, we shall also see that the strength of this repulsion is 
exactly the one that gives the correct equation of state ( p.07| ), (|49|) at small densities: 

ub = 32NfF' (^;, - j + . . . . (108) 

This is a very non-trivial check of the consistency of our effective Lagrangian! 

Since we are working at very small densities, we can expand the Lagrangian (^9|) up to 
next (quartic order) in the Goldstone fields. The quartic terms will describe the interaction 
energy of the Goldstones due to two-body scattering, 

^int = -^Tr [n, d,U] [n, d^U] - ^Trn^ + . . . . (109) 

We are now working in Minkowski space and the extra minus sign in front of the first term, 
compared to (|29|) , is a consequence of this. The presence of the time derivatives should not 
confuse the reader. Since we are dealing with a non-relativistic system, /i — <^ m^, 

we only need the leading term in the time dependence of the relativistic fields, i.e. doll —>■ 
izmTrll. The spatial derivatives give only a subleading contribution. The dominant effect is 
that of the s-wave amplitude. 

We now consider the ground state of such a system with the baryon density fixed (unlike 
before, when we fixed fi). Most of the particles will occupy the lowest energy, i.e., zero- 
momentum, level (the population of the excited levels is the next order effect). These 
particles form a condensate, whose amplitude we denote by 11. The energy of the system 
is different from zero only because of the interaction, and, as a function of 11, it is given 
to us by ( |109| ). Looking at the block decomposition of 11 (p!6| ) we note that only Q and 
(diquark and antidiquark) fields carry non-zero baryon number. Therefore, we set P = in 
our case. Thus, from ( p.09| ), we find for the dependence of the energy on the magnitude of 
the Q condensate 

2 

S,^, = ^Tr(gQt)2 + . . . ^ (110) 

where we neglected subleading non-relativistic corrections. 

Now we need to calculate as a function of Q. Since Q carries baryon charge 2, we can 
writeQ 

nB = 2 TiiQ^doQ - Qid^Q^) = Am^TiQQ^ + ..., (Ill) 

where we neglected subleading corrections. Since all flavors of diquarks are equivalents^ the 
matrix Q is proportional to / for /3 = 1, or to 1 for /? = 4. Thus we can relate traces in 

~ 

2F2TrBI]t5oE = 2iTiBUdoll, which coincides with (pTT|) 



Another way to see (till) is to use the fact that ub — —dC/d^i. Thus, from (WM we find that ns = 



There is only one flavor of diquarks for (3 — 1, Nf — 2. For larger Nj the flavor SU{Nf) symmetry 
ensures the equivalence of all diquarks. 



33 



(110) and (|1 1 1|) by Tt{QQ^)^ = (TiQQ^)'^ /Nf. Therefore, we find for the dependence on 
of the vacuum energy at fixed ub, 

^ ^B ^ Ub n\ 

^vac = + ^int = -^m^ + g^^,^^ + . . . . (112) 

The first term is the rest energy of the diquarks. This determines the equation of state 

The fact that the equations of state ( |108|) and (|113|) are identical is a very nontrivial property 
of the effective Lagrangian (0). It is intimately related to the symmetries (global and local) 
which determine the form of this Lagrangian. 

This calculation also shows that the diquarks Bose-condense and thus form a super- 
fiuid. Since the diquarks are charged we can call this phase a superconducting phase. 



14 Conclusions and discussion 

In this paper we have derived the low-energy effective Lagrangian for certain QCD-like 
theories at finite chemical potential fi. It contains the same number of phenomenological 
parameters as at = 0: the pion decay constant F and the vacuum value of the chiral 
condensate {ipip)o. Using this effective Lagrangian we have completely determined its low- 
energy properties. In particular, we have established the dependence of the ground state, the 
condensates, and the masses of the excitations on the chemical potential /i, the bare quark 
mass m and the diquark source j. The theories to which our analysis applies include: (i) two- 
color QCD with quarks in fundamental color representation, and (ii) any-color QCD with 
quarks in adjoint color representation. The unifying feature of these theories is that the low- 
energy excitations, i.e. the Goldstone bosons of spontaneously broken symmetries, include 
baryons in the form of diquark states. Since the chiral Lagrangian describes the dynamics of 
such low-energy excitations, we can use the effective theory to describe phenomena associated 
with the condensation of diquarks. 

Since we are working only to the lowest order in Chiral Perturbation Theory, we must 
remember that our approach has a limited region of applicability. First of all, the momenta 
of the particles, described by the Lagrangian ( ^2]) must be smaller than the scale M of the 
masses of the non-Goldstone excitations, such as vector mesons, for example. In other words, 
the expansion parameter of the chiral perturbation theory is p/M. The bare quark mass 
is also included to leading order in m. It is obvious from (^) that m ~ ~ p"^ in the 
momentum power counting. The chemical potential fi is an external parameter with the 
dimension of mass. The interesting behavior of our theory is at values of /x ~ ~ p, and 
fi is counted as order p in the momentum power counting. Our effective theory thus applies 
as long as and fi are much smaller than M. 

Our effective theory predicts a phase transition at a value of fi = It is important 

to note that this transition occurs within the domain of applicability of the effective theory, 
as long as is small, as discussed above, unlike, for example, the transition in QCD with 
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three colors. The transition at /i = is a robust prediction of our effective theory. The 

existence of such a transition is confirmed by lattice Monte Carlo simulations ^. 

The fact that our effective Lagrangian ( ^21) predicts the transition at a value of /i equal 
to is a direct consequence of the symmetries of the theory. The coefficient of the 

/^-dependent term in the effective Lagrangian, responsible for the transition, is not an inde- 
pendent phenomenological parameter, but is fixed rigidly by the local symmetry ||^. This 
symmetry corresponds to the conservation of the baryon charge (among other generators of 
SU{2Nf)). This ensures that the charge of a composite state in the effective theory is an 
additive sum of the charges of the constituents (non-renormalization). Indeed, the factor 
two in the relation /iq = is the baryon charge of the diquark. This is in agreement 

with the value of /i, below which no transition can occur at zero temperature, given by the 
minimum value of the mass per baryon number among all baryons in the theory (see, e.g., 
||3^). In three-color QCD this would correspond to a value of /iq approximately equal to 1/3 
of the nucleon mass (offset by the binding energy of nuclear matter) .[^ 

Another way of looking at the /i-dependence in the theory is the following. In the 
microscopic theory, given by the Lagrangian (PB|), fi enters as (a timelike component of) an 
Abelian gauge potential. As such, it can be completely removed from the Lagrangian by the 
time-dependent gauge transformation 

ijj e^^^il) and i)* e'^^*, (114) 

where r is the Euclidean timej^. However, this does not mean that the partition function 
does not depend on /i. Let us consider first the case of finite temperature. In this case, the 
boundary conditions in the Euclidean time direction for the quarks change after (|114|) : they 



are no longer antiperiodic. The dependence on in the partition function comes entirely 
from the boundary conditions on the fermion fields 

Mr=i/T = -e^/^^|.=o, and r\r=i/T = -e-^'^r\r=o- (115) 

The /i-independent Lagrangian with such boundary conditions is completely equivalent to the 
Lagrangian ( P^ and usual antiperiodic boundary conditions. This fact can be conveniently 
used in lattice simulations. The corresponding effective theory can also be defined with 
all /i-dependence in the boundary conditions. The boundary conditions for a given effective 
(composite) bosonic field should then read 0(l/r) = exp(6/i/T)0(O), where h is the baryon 



charge of the field 0]. In the case of the matrix valued field S we have 



T.l=i/T = exp {^^B^ S|,=oexp {^^B^^ , (116) 

where B is the baryon charge matrix (^). 

In the limit of zero temperature, our intuition suggests that the dependence of the par- 
tition function on boundary conditions should weaken and disappear. This is, however, not 
completely true for the boundary conditions such as ( |115|) due to their singular nature in 
the limit 1/T — >■ oo. There is indeed an interval of /i: — /io < /i < /xq, where the partition 

i2\Ye remind the reader that we measure the baryon charge in units of the U{1)b charge of a single quark. 
^•^ Note that, in Euchdean formulation, the ij} and ij}* are independent variables and the global flavor 
symmetry group is in fact Gl{2Nf) p^ . 
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function does not depend on /x. However, outside this interval, the singular nature of the 
boundary conditions starts playing a role: a phase transition occurs and a //-dependence 
appears. 

Why is the approach of this paper not directly applicable to real three-color QCD? The 
main reason is that at /i = the effective theory described by the chiral Lagrangian does 
not contain excitations with non-zero baryon number (pions do not carry baryon charge). 
It is easy to see that, as a consequence, applying either the method of local symmetry of 



Section |6]^ or the gauge transformation described in this section, one finds no dependence 
on /i in the effective chiral Lagrangian of QCD. Related to that is the fact that the value of 
/io ~ 'Tinucieon/S is large, and the approach of this paper, based on an expansion in small 
will not reach it. 

As we emphasized, the theories we considered are related by the fact that the fermion 
representations are pseudo-real. This is a consequence of the antiunitary symmetries of 
the Dirac operator (0), (^). We identified two such symmetry classes, distinguished by the 
Dyson index j3 = 1 and j3 = 4. The effective theories for these two classes of theories 
are very similar and dual, or complementary, to each other, from the point of view of the 
residual global symmetries. In the diquark (finite density) phase, for the (3 = 1 theories 
the residual flavor symmetry is given by Sp{Nf), while for the (3 = 4 case it is SO{Nf). 
The excitations form multiplets which correspond to symmetric or antisymmetric second 
rank tensor representations of these groups (see Figs and Tables |l]Jl). One can see 
that the /? = 1 and (3 = 4 cases "mirror" each other with respect to Sp{Nf) ^ SO{Nf) 
and symmetric ^ antisymmetric. On the other hand, the Dirac operator for three-color 
QCD with fundamental quarks does not have any antiunitary symmetries, and the fermion 
representations are complex. This case falls into the third remaining Dyson class with the 
index P = 2. 

It would be interesting to apply the approach of this paper to lattice theories with pseu- 
doreal fermions, for example to two-color QCD with fundamental quarks. The symmetries 
of such theories are different from their continuum counterparts and were analyzed in 0. 
In particular, the transition to continuum limit may turn our to be nontrivial. This can be 
related to the antiunitary symmetry of the lattice Dirac operator for staggered fermions in 
fundamental representation 

where rix^fj, = (_ i)^'i+'--+^m-i^ and U are SU{2) color matrices. The antiunitary symmetry of 
this lattice Dirac operator is given by T2D = V*T2, or \I),T2K] = 0. Since {T2KY = —1, we 
conclude that such Dirac operator belongs to the class (3 = 4. However, in the continuum 
limit, the Dirac operator must be in the class (3=1. Such an observation was also made 
in ^ |2^, from the point of view of global symmetries and their breaking, where it 



was pointed out that it is not yet known how the apparent pattern of SU{2Nf) 0{2Nf) 
breaking becomes SU{2Nf) —>■ Sp{2Nf) in the continuum. Even though the symmetry 
may modify many of the details of our analysis, when it is applied to lattice theories, some 
features should be robust. Such would include the phase transition at /i = 171^^/2, the relation 
similar to ('?/''?/') ^ + (V'V')^ = const between the chiral and the diquark condensates, the linear 
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dependence of on near the transition, the existence of several branches in the spectrum, 
similar to Figures |l|J|,|, and many other qualitative features. 
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Appendix 

In this appendix we analyze the representations of the remaining symmetry groups shown 
in Table |I]for /5 = 1. They are all different subgroups of Sp{2Nf). 

Under the symmetry group Sp{2Nf) the Goldstone fields S transform as 

s ^ vj:v^. (118) 

With S parameterized as U^U^ (with U G SU{2Nf)) and KScK^ = (the 2Nf x 
2Nf antisymmetric unit matrix is denoted by Sc) it follows that the generators transform 
according to 

XS, ^ VXV-'K = VXT.,V^ or (XS,),, ^ V,uV,i{XT.,)u. (119) 

From the transposition relation ([T^ ) it follows that (XSc)"^ = — XSc. If the symmetry 
group is Sp{2Nf) the generators transform according to an antisymmetric rank two repre- 
sentation of Sp{2Nf ). The degeneracy is thus given by 2Nj — Nf — 1. If the symmetry group 
is SU{Nf) X U{1) the symmetry transformation is given by 

(120) 

the P-type generators transform according to UiP'^Ui^ and the Q-type generators as 

Q — > UiQUf . Since the Q are antisymmetric they transform according to an antisymmetric 
rank two representation of SU{Nf). The fields transform according to the conjugate 
representation. The dimension of both representations is equal to Nf{Nf — l)/2. Since the 
P-type generators are traceless, the degeneracy is given by Nj — 1. 

In the diquark condensation phase the symmetry group Sp{2Nf) is with respect to the 
rotated antisymmetric unit matrix, S^^, with symplectic transformations defined by 

VE^V^ = (121) 




37 



According to the argument at the beginning of this appendix, the transformation properties 
of the generators are given by 



{XaT,a)ij — > VikVji{XaT,a)kl, (122) 

where the Xa are the rotated generators defined by 

X„ = VaXV-\ (123) 



and Va = exp 10X2/2 defined in 

With Sp{Nf) as symmetry group the symmetry transformations are given by (|120|) with 
U* = —lUiI (the Nf X Nf antisymmetric unit matrix is denoted by /). Since in this case 
{V, X2} = and using that S = VaUT^cW^V^ (see eq. (|6^)) we find that the representations 
can be discussed in terms of the fields UT^JJ'^ with the famihar PQ-block structure of the 
generators. The P-type generators thus transform as I UiP^ lUi . The symmetric and 
and antisymmetric components of PI corresponding to PsI and PaI ■, respectively, transform 
independently. The dimensions of the representations are given by Nf{Nf+l)/2 and Nf{Nf — 
1) — 1, respectively. The generators Q and IQ^I transform in the same way, and thus 
the linear combinations Q and transform in the same way. Since Q is antisymmetric 
they transform according to an irreducible rank two representation with dimension equal to 
Nf{Nf - l)/2. 

Finally, for symmetry group Sp{Nf) x Sp{Nf) the symmetry transformation is given by 



v= r,; ' , (124) 




with both Ui and U2 symplectic transformations. In this case it is imperative to consider the 
rotated generators. However, we only need the rotated generators for /i ^ 00. In terms of 
the block structure of the unrotated generators, we obtain after a straightforward calculation 
for a = 7r/2, 

_(-IPaI + iQiI iIPs + Qr\ 

Therefore the combination ilPs + Qr transforms as 

iIPs + Qr^ U^IPs + Qr)U^. (126) 

Indeed, this transformation corresponds to the Young tableaux given in upper right corner 
of Fig. 1. Notice that for /i — >■ cxd we have that Qr. The 11- and 22-blocks of the 

matrix ( |125| ) correspond to two different linear combinations of Pa and Qj. One combination 
transforms as a rank two tensor with respect to f/i, and the other combination as a rank 
two tensor with respect to U2- Since both combinations are antisymmetric with respect 
to transposition, they transform according to the Young tableaux given in the lower right 
corner of Fig. 1. The additional singlet terms arise because the irreducible representations 
are traceless. 

In case of Sp{Nf) symmetry, U2 = = —lUiI, the transformation properties of P5, 
Pa, Qr and Qj can be obtained from ( |125| ) by combining the transformation properties of 
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the diagonal blocks and of the off-diagonal blocks. The results are in agreement with the 
discussion in the paragraph following eq. (|123|) . 

The same analysis can be performed for (3 = 4. In fact, because the matrix / is absent, 
this case is somewhat simpler, and we leave it as an exercise to the reader. 
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